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I. INTRODUCTION 

Supercooled liquids and glasses have very large re- 
laxation times and complex relaxation mechanisms, but 
their structures appear disordered and unremarkable [T]- 
3 . This combination is surprising and rather mysteri- 
ous, especially because several recent studies [1HB] indi- 
cate that if a relaxation time increases sufficiently rapidly 
on cooling then this must be accompanied by the devel- 
opment of some kind of structural order. The idea of 
growing amorphous order is that the diversity of amor- 
phous states, as quantified by the configurational part of 
the entropy, decreases at low temperatures. This leads 
to increasing static correlation lengthscales, which can be 
measured using point-to-set correlations [41410J . If these 
length scales are large, the system may be localized into a 
single amorphous state by fixing the positions of a small 
fraction of the particles [TTJ [T^] . This localization occurs 
at temperatures above the glass transition temperature, 
so the pinning procedure bypasses the challenging task 
of thcrmalizing systems at very low temperatures. 

Here, we investigate amorphous order in two finite di- 
mensional interacting spin systems. These are the square 
plaquette model (SPM) [13] and the triangular plaque- 
tte model (TPM) [TJ]. They both have growing relax- 
ation times at low temperatures, but two-point thermo- 
dynamic correlation functions do not indicate the pres- 
ence of any growing length scale [15], and the models 
do not have phase transitions at any finite temperature. 
Additionally, structural relaxation at low temperatures 
exhibits strong dynamical heterogeneity and growing dy- 
namic lengthscales [TBI HZ] • These features mimic those 
of supercooled liquids, but the models are numerically 
and analytically more tractable than off-lattice particle 
systems. This makes them useful models for studying 
generic features of glassy systems [TfjHT5] , 

Several methods have been proposed for characteriz- 
ing amorphous order. A prominent recent example is 
to fix particles everywhere except within a small cavity, 
and then to study the motion of particles within this 
cavity [U [7] [5] US [2D]- Other possible measurements 



of amorphous order involve freezing spins in a different 
geometrical arrangement [9[ 111] , or direct measurement 
of the diversity of amorphous packings [H HU] ■ Here we 
follow [TTJ [TH |2"TH2"3"] and fix a fraction / of particles 
(or spins), distributed randomly through the system in 
the positions (or orientations) that are representative of 
thermal equilibrium. In supercooled liquids, a thermo- 
dynamic phase transition on varying / was recently pre- 
dicted [12j . based on an analysis within random first- 
order transition (RFOT) theory [2H [25]. Another re- 
cent analysis based on mode-coupling theory predicts a 
dynamic singularity on varying / |26j . Numerical stud- 
ies [TT] |2"TH2"3"] have mainly addressed the behaviour at 
low to moderate pinning fraction /, leaving open the 
question of the existence of phase transitions at larger 

Here, we explore the effect of random pinning in the 
plaquette models, concentrating on the behaviour at the 
relatively large pinning fractions where phase transitions 
are predicted to occur in supercooled liquids. We show 
that the plaquette models do not exhibit thermodynamic 
phase transitions on varying / at any finite temperature. 
However, we do find a well-defined crossover from 'bulk' 
behaviour (where the system explores a large number of 
configurations), to 'frozen' behaviour (where the system 
remains trapped in a region of configuration space) oc- 
curring at some finite fraction /*. As the temperature 
is reduced, the crossover occurs at an increasingly small 
value of /*, and becomes increasingly sharp. 

From previous studies of point-to-set correlations in 
closed cavities for plaquette models [T7] , it is known that 
these models do have growing amorphous order at low 
temperatures, but their behaviour is not consistent with 
RFOT. The results we present here therefore illustrate 
a theoretical scenario, alternative to mode-coupling and 
RFOT treatments, for understanding the behaviour of 
model supercooled liquids with randomly pinned parti- 
cles. 

After defining our models in Sec. |Tt[ w e show numer- 
ical results for this crossover in Sec. Mil In Sec. IIVI we 
investigate the length scales that characterize amorphous 
order in the system, including their scaling with temper- 
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ature. In Sec. [V]we use analytic arguments in the SPM 
to study the small-/ behaviour, and the behaviour near 
the crossover. Finally in Sec. |VI| we discuss the interpre- 
tation of our results, including their relation to previous 
theoretical analysis. 



II. PLAQUETTE MODELS 
A. Definitions 

The plaquette models consist of Ising spins, Sj = ±1, 
on a lattice, interacting through either three- or four- 
body terms. They evolve dynamically through single 
spin-flips with Metropolis rates, Wi = min(l, e _/3AEi ), 
where AEi is the change in energy of the system on flip- 
ping spin i and ft = 1/T is the inverse temperature. 

The energy of the square plaquette model (SPM) is 

E = -- ^2s 1 s 2 s 3 s 4 , (1) 
□ 

where the sum runs over plaquettes of the square lat- 
tice [131 US]- It is useful to define c = c~^, since length 
and time scales typically show power law scaling with 
c in the low temperature limit (c — > 0). The model 
has a dual representation in terms of 'defect' variables, 
rip = (1 — S1S2S3S4V2, which are associated with the 
plaquettes p of the square lattice. In terms of these 
defect variables, the thermodynamic properties of the 
system are those of a non-interacting lattice gas with 
(n p ) — (1 + C 9 ) -1 approximately equal to c at low tem- 
peratures. We use square systems with periodic bound- 
aries. 

The energy of the triangular plaquette model (TPM) 

is 

£ = ~5>is 2 s 3 , ( 2 ) 
v 

where the sum runs over downward pointing plaque- 
ttes of the hexagonal lattice [2]. In this model, the 
defect variables are associated with these plaquettes, 
n p = (1 — siS2S 3 )/2. We use rhombus-shaped systems 
with periodic boundaries. We always take the linear sys- 
tem size to be a power of 2, which minimizes finite size 
effects. 

In the defect representation, both the SPM and the 
TPM have strong similarities with kinetically constrained 
models [27] . In both systems, the thermodynamic free 
energy per spin (in the limit of large systems) is simply 
— Tlog[2cosh(l/2T)]. This is a smooth and unremark- 
able function of temperature down to T = 0. However, 
the dynamical properties of the models, to be described 
shortly, are not trivial. This combination of simple ther- 
modynamic properties and complex dynamics is consis- 
tent with the dynamic facilitation picture of the glass 
transition [28] . 



B. Square plaquette model at / = 

Before considering the effect of pinning a fraction of 
spins, we first review some relevant results for the SPM 
with / = 0, see Refs. [TH HE] f° r more details. Sym- 
metry under global spin reversal implies that (s,) = 0, 
and all non-trivial two-point and three-point correlation 
functions also vanish. The simplest non-trivial correla- 
tion functions involve four spins. For spins a, 6, c, d lying 
on the four vertices of an x x y rectangle then 

g a (x, y) = (s a s b s c s d ) = tanhfj3 /2) xy . (3) 

All other non-trivial four-point correlations vanish. 

The function g s (x, y) has scaling behaviour at low tem- 
peratures. Full details are given in Sec. |IVB] below, but 
the essential points are that correlations are anisotropic 
and the model supports two correlation lengths at low 
temperatures. The shorter length is £0 ~ c -1 / 2 , and 
there is also a longer length £ ~ c" 1 associated with the 
correlations along the lattice axis. We note that previous 
analysis of static correlations considered only the shorter 
of these two correlation lengths [16] . 

Point-to-set correlations were considered in Rcf. [17j . 
for the case of an isotropic (square) cavity of linear size 
r. For r < £0 ~ c^ 1 / 2 , the configurational entropy of the 
cavity is close to zero. For £0 < r % £ ~ c _1 , the cavity 
has strong finite size effects but its configurational en- 
tropy is non-zero and the spin-spin autocorrelation func- 
tion decays to zero at the centre of the cavity. This im- 
plies that these larger cavities are no longer frozen in a 
single amorphous state, so we identify £ as the point-to- 
set length [17] , 

Dynamical observables in the square plaquette model 
also show scaling behaviour. The relaxation time for 
spins is r w c~ 3 ~ e 3 / T , consistent with the behaviour 
of 'strong' glass- formers [TS] . Energy-energy correlations 
decay have relaxation times longer than t at low tem- 
peratures, although they also have Arrhenius scaling. 
Four-point dynamic correlation functions, also discussed 
in more detail below, are anisotropic as well, with cor- 
relations being strongly localized along the axis of the 
square lattice, and having spatial on-axis extension of 
range £ ~ C 9 = c _1 16]. 

C. Triangular plaquette model at / = 

Static correlations in the TPM were discussed in 
Ref. [13] where it was shown that (s,) = and (stSj) — 
5ij, as in the SPM. The relaxation time in the TPM 
has 'fragile' super- arrhenius scaling r ~ e 1 l^ T log3 ) [IT] . 
which arises from a hierarchy of mechanisms whose char- 
acteristic timescales increase as the logarithm of their 
associated lengthscales (see also [2"M3"T] ). Static corre- 
lation lengths, point-to-set correlations, and dynamical 
four-point correlations are considered in Ref. |17j . where 
it was found that the scaling of these functions all de- 
pend on a unique correlation lengthscale, £ ~ c _1 / df , 
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where d{ = log 3/ log 2 ss 1.585 is the fractal dimension 
of the Sierpinski triangle. It was also found recently [32] 
that the TPM supports unusual phase transitions if it is 
constrained to have a fixed non-zero magnetisation, but 
we do not discuss that case here. 



III. EFFECT OF RANDOM PINNING 

We now turn to our main results, which correspond 
to the following thought-experiment [TT | [T2 | I2T1 122) 126 ] 
155] . We select a reference configuration from a thermally 
equilibrated system and we instantaneously 'freeze' (or 
'pin') at time t = the state of a finite set of spins. 
Each spin is frozen with probability /, independently of 
all other spins. The spins that are not frozen evolve with 
Monte-Carlo (MC) dynamics for some time t, after which 
various measurements are performed. 



A. Correlation functions and susceptibilities 

To describe these measurements, we first establish our 
notations and define the correlation functions and suscep- 
tibilities that we will consider. The system consists of V 
spins Si = ±1, with i = 1, . . . , V. To describe whether a 
spin is frozen, we introduce the binary variables /j, with 
fi = 1 if spin i is frozen, zero otherwise. We use angle 
brackets (• • • ) to denote averages which run over the ref- 
erence configuration, the choice of frozen spins, and the 
MC dynamics. We also use (• • • )o to denote a 'bulk' ther- 
mal average, i.e. in the absence of any pinned particles 
if <»• ' 

It is convenient to define the autocorrelation of the 
mobile spins. For spin i, this is defined as ai(t) = 
(1 — /j)si(t)si(0), and it takes values 0, ±1. The num- 
ber of mobile (unfrozen) spins is N m — 5^(1 — /i), which 
is a fluctuating quantity within our analysis. We also 
define the (extensive) overlap Q(t) = J2i a i(t) anc ^ the 
(intensive) autocorrelation function 



C(t) = 



(Q(t)) _ (ai(t)) 



1-/ 



(4) 



The four-point susceptibility quantifies the strength of 
spontaneous fluctuations of the overlap, 



X4(t) - 



(SQ(t) 2 ) = V,/y.u 

(N m ) l-f ' 



(5) 



and it is related to the volume integral of the four-point 
correlation function defined as 



94,ij = (Sai(t)8aj(t)). 



(6) 



We use the notation SO = O — (O) for thermal fluctua- 
tions throughout this article. 

In addition to the usual four-point correlations, we in- 
troduce three-point correlations and susceptibilities. We 



measure how the positions of the pinned spins are corre- 
lated with the states of mobile spins through the corre- 
lation 



(7) 



The ratio (A^//) measures the change in the average of 
di if spin j is assumed to be frozen. Hence, 



dC{t) 



1 



/(i-/) a £T 



(8) 



e^) _1 , so that 
follows. 



To prove this relation, one writes / = 
(<»*(<) E^/i) = -|>i(*)>, and Eq ^ 

It is interesting to note that X# (t) Is" an example of a 
three-point correlation function, and that dC g^p is there- 
fore a three-point susceptibility, in the spirit of the dy- 
namic functions introduced in Refs. [5H |55] . Following 
Refs. [34-36] we can obtain a lower bound A%4 on the 
four-point susceptibility, Xi{t) > ^X&{t)i which involves 
the three-point susceptibility defined above as 



A X4 (<) = / 



This relation is proved in Appendix|A]where we also show 
that a sufficient condition for saturation of the bound 
[i.e., Xi(t) — AX4(*)] is that for a given choice of frozen 
spins, the autocorrelations of the mobile spins are inde- 
pendent of each other and depend separately on the indi- 
vidual fj . When this condition is obeyed, then the four- 
point correlation function may be expressed as a convo- 
lution of Xij with itself: 



/( I^7)E^W^W- 



(10) 



We refer to Appendix [A] for more details on this partic- 
ular point. 

Of particular interest is the limit of large time, t — > oo, 
in which correlation functions become independent of the 
dynamical evolution of the system and may be calcu- 
lated by equilibrium statistical mechanics. In this limit, 
we drop the time argument on our correlation func- 
tions, writing, as t — > oo: C(t) — > q, g<i,ij(t) — > 94,,ij, 
Xij(t) — > X^, and Xi(t) ~^ Xi- These limiting quan- 
tities are given by 'static' correlation functions. A gen- 
eral recipe for calculating them is given in Appendix [A] 
which follows a similar analysis for off-lattice particle sys- 
tems |37) . For small /, in particular, they may be calcu- 
lated in a series expansion, see also Sec. [V] 

B. Numerical results 

An overview of the influence of pinned spins on the 
SPM is shown in Fig. [T] In Fig. [IJi we show that as / is 
increased from zero, the long-time limit of C (t) increases 
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FIG. 1: Overview of behaviour of the SPM with pinned spins, (a) Time-dependent correlation function C(t) at inverse 
temperature /3 = 3. (b) The long-time overlap q — lmit-joo C(t) as a function of pinning fraction / for four temperatures (the 
symbols are the same as those in panel e). (c) Time-dependent four-point susceptibility Xt{t) at /? = 3. (d) Long-limit limit 
of X4, for various temperatures, (e) Behaviour of the relaxation time on varying / at /3 = 3. (f) Behaviour of the maximal 
relaxation time r max = max/ t(/), normalised by the bulk relaxation time to = r(f — 0), as a function of inverse temperature. 
In practice r max « r(/*). The straight line is r max /ro oc e 13 which would imply r max oc e 4/3 . It is clear from the upwards 
curvature of the numerical data that r max has a super-Arrhenius temperature dependence. 



monotonically, since the frozen spins tend to maintain the 
system close to its initial state, preventing full relaxation 
of the autocorrelation function of the mobile spins. Thus, 
we confirm that by randomly freezing spins, the system 
crosses over from 'bulk' ergodic relaxation for / = to 
a nearly 'frozen' state at large /. Our goal is now to 
characterize this crossover further. 

In Fig. [T]d, we show the evolution of q — q(f,T) with 
the fraction of pinned spins and temperature. As ex- 
pected from Fig. [T^,, the static value of the overlap in- 
creases monotonically with / at any given temperature. 
More interesting is the temperature dependence of the 
q(f, T) curves. We find that q increases rapidly with de- 
creasing temperature at constant /. This implies that for 
lower temperatures, a smaller amount of random pinning 
is required to localize the system in a single state. The 
interpretation is that the system has a greater degree of 
amorphous order at low temperatures. 

Looking more closely at the /-dependence of q, the 
data in Fig. [T|d indicate that q(f,T) has an inflexion 
point at a characteristic value of the pinning fraction, 
/ = /*, so that the susceptibility (dq/df) is small both 
for small / and for large /, with a well-defined maxi- 
mum at f*(T). Thus we find that the bulk-to- localized 
crossover obtained with random pinning can be located 
by measuring the derivative of the static overlap. How- 



ever, anticipating the discussion in Sec. |VI[ we note that 
(/(/, T) is a smooth function of / with no sign of the 
sharp discontinuity that would be observed at a first or- 
der phase transition 

We next turn to fluctuations of the overlap, which 
we quantify via the four-point susceptibility xa{^)- I R 
Fig. [TJ;, we show the time evolution of Xi(t) f° r differ- 
ent values of /, at constant temperature. For / = 0, 
the susceptibility X4,(*) uas a peak for t s» r, as usual in 
glassy systems 10J . However, two features emerge when 
/ is increased. First, the time dependence changes dra- 
matically: the maximum in Xiif) shifts to longer times 
until, for large /, the susceptibility Xi{t) is monotoni- 
cally increasing and saturates to a plateau at long times. 
This long time limit corresponds to the static suscepti- 
bility xa — XiifiT), proportional to the variance in the 
overlap q(f, T). This increasing static susceptibility indi- 
cates that deviations between the final configuration and 
the initial (reference) configuration appear by coopera- 
tive processes involving many spins. Just like (dq/df), 
the static susceptibility Xi goes through a maximum at 
the characteristic pinning fraction /*. 

Physically, the interpretation of the behaviour of Xi{t) 
is as follows. When / <C /*, the system is in the 
bulk regime and easily escapes from the reference state 
through a process that is not very different from bulk 
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relaxation. In this case, Xi(t) is large near t « r, but 
it is small at long times since initial and final states are 
very different. When / 3> /*, by contrast, there are so 
many frozen spins that the system is very constrained 
and few spins can relax. While the overlap is large, its 
fluctuations are necessarily quite small. For / ~ /*, the 
number of frozen spins is just large enough to maintain 
the system near its initial state, and the overlap exhibits 
stronger fluctuations because the system 'hesitates' be- 
tween both possibilities ("should I stay or should / go"). 
We discuss the spatial structure of these correlations in 
Sec. HV] below. 

We show in Fig. [TJi the evolution of the static suscep- 
tibility Xi with / for different temperatures. As antici- 
pated, the susceptibility goes through a maximum whose 
location and amplitude are strong functions of the tem- 
perature. On going to lower temperature, the peak of \i 
remains located near / = /* so it shifts towards smaller 
values of /; the amplitude of the peak increases rapidly, 
and its width decreases. Thus, the crossover between 
bulk and localized behaviours becomes sharper and more 
pronounced at low temperatures. A relevant conclusion 
for supercooled liquids is that the data in Fig. [TJl show an 
increasing static susceptibility that measures the growth 
of amorphous order, but these data are obtained without 
any a priori knowledge of the many-body correlations 
that are responsible for this order. Thus, the random 
pinning procedure is a generic way to measure a grow- 
ing static susceptibility in liquids approaching the glass 
transition, and offers a thermodynamic alternative to the 
measurement of relevant length scales via dynamic het- 
erogeneity. 

Turning to the dynamic behaviour in the presence of 
random pinning, we define a relaxation time r = r(/, T) 
from the time decay of C(t), via c ^~ g = 1/e. We com- 
pute r using the data shown in Fig. [lfi and show the re- 
sults in Fig.[T^. We also find that r has a non-monotonic 
behaviour, the relaxation being slowest near /*. While 
the maximum is not very pronounced in Fig. [l^, we 
show in Fig. [If that the ratio r(/*,T)/r(/ = 0,T) in- 
creases when temperature is reduced. This is consis- 
tent with the presence of increasingly cooperative re- 
laxation mechanisms in the presence of pinned spins. 
While the SPM without pinning has 'strong' glass scal- 
ing, t(/ = 0) ~ e 3 / T , we find that r(/*) increases in a 
super- Arrhenius (fragile) fashion. This indicates that the 
relaxation mechanism for the system near /* is different 
from the bulk mechanism at / = 0, presumably because 
the frozen spins act to frustrate relaxation of the mobile 
ones. 

We expect many of the results shown for the SPM to be 
quite generic in glassy systems. Certainly, the behaviour 
in the TPM is similar. Figure [2] shows data for C(t) 
and %4(t) in the TPM at the representative temperature 
(3 — 3. As before, C(t) exhibits a plateau at long times 
which increases with /, while the static susceptibility %4 
is maximal at some / = /*. Consistent with the SPM, 
we also find that /* decreases at low temperature, that 





'$ = 3 


/= o.oo -H- v 




0.05 -*- 




0.10 




0.15 -— 




0.20 ■ 





10° 10 1 10 2 10 3 10 4 10 5 

t 




10° 10 1 10* 10 3 10 4 10 5 



FIG. 2: Effects of pinning in the TPM at = 3. (a) Time- 
dependent correlation function, (b) Time-dependent four- 
point susceptibility. 

the maximum of the static susceptibilities increase, and 
that the time scales at /* increase. These results resem- 
ble strongly the ones shown in Fig. [T] and are therefore 
not shown, for brevity. However, the /-dependence of 
the relaxation time r is weaker in the TPM than in the 
SPM: the ratio r max /ro does increase systematically on 
decreasing temperature but it takes values in the range 
1 — 4 while To varies over nearly four orders of magnitude. 

We emphasize that all of the results presented in this 
article are obtained in large systems and we have checked 
that they are free from finite size effects. In particular, 
susceptibilities and relaxation times have maxima at /* 
but the maximum values remain finite even when the 
thermodynamic limit is taken, V — > oo. Thus, there are 
no diverging correlation times or correlation lengths in 
this system for any finite / or T, nor is there any phase 
transition. However, we find that correlation times and 
susceptibilies have sharp maxima along a line /* (T) in 
the (/, T) phase diagram. 



IV. SCALING OF LENGTHS 

We have shown that varying / in the SPM and TPM 
reveals crossovers at /*, associated with maxima in sus- 
ceptibilities and in relaxation times. We now discuss how 
these features can be related to correlation lengths in 
these systems. In particular, we focus on the scaling of 
these lengthscales at low temperatures. 
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FIG. 3: Snapshots of the propensities p t in the long-time 
limit, (a-c) SPM at /3 = 3, varying /. At this temperature 
/* « 0.18. (d-f) TPM at /3 = 3 for which /* « 0.15. Pinned 
spins are black while unpinned spins are color-coded with dark 
blue for pi w 1, pale blue for pi w and green for pi ~ —1. 
On increasing /, the spins become polarised so that p; w 1 in 
most cases. The correlations associated with the propensity 
are strongest near /*, consistent with \i being maximal. 



A. Visualisation of spatial correlations 

It is instructive to visualise the spatial fluctuations 
that appear as a result of the random pinning. To this 
end, we consider the dynamic propensity [351 [35]. (Com- 
pared to visualising the autocorrelations at directly, the 
propensity provides continuous functions rather than bi- 
nary ones, and this yields images that better differenti- 
ate between regions where relaxation is frustrated by the 
frozen spins and those where relaxation can occur.) 

To calculate the propensity, we take a single repre- 
sentative reference configuration, s A , in which a specific 
set of spins are frozen, and we run several long MC tra- 
jectories starting from it. We calculate the autocorre- 
lation a,i(t) = Si(t)si(0) for each unfrozen spin in each 
trajectory and we average over the trajectories to ob- 
tain the (site-dependent) propensities pi(t) = di(t). For 
large times these propensities approach limiting values, 
Pi = Pi(t which depend on the reference config- 

uration s A but not on the time t. The propensities are 
therefore static on-site quantities characterizing the de- 
gree of freezing of spin i for a given realization of the 
random pinning and a given reference configuration. 

Representative results are shown in Fig. [3J where 
pinned spins are shown in black, and a blue-green colour 
coding describes the propensity. Sites for which the 
pinned spins cause the configuration to remain near its 
initial state have pi ps 1 (dark blue) while those where the 
pinned spins have little effect have pi — (light blue) . If 
the frozen spins cause Sj to become polarised in the op- 



posite direction to sf~, one finds pt < (green). 

At / = then pi = for all i and the system is ho- 
mogeneous. As / starts to increase (left panel), then 
the system acquires regions where the spins become po- 
larised and cannot relax any more (coloured dark blue 
in Fig. [3]) . This seems to occur in small, isolated regions 
whose size increases with /. For / ~ /* (middle panel), 
these polarized regions percolate throughout the system 
and spatial fluctuations of the propensity occur over a 
large lengthscale. This yields snapshots where large re- 
gions are strongly polarised while others are unaffected 
by the pinning. Finally, as / increases further above /* 
(right panel), most spins are strongly pinned, and only 
few small regions exist where motion remains possible. 
In this regime, the system is strongly localized near the 
reference configuration. 

An important observation is that the lengthscales asso- 
ciated with these correlations observed near /* are much 
longer than the typical distance f" 1 ^ 2 between pinned 
spins. This is apparent in Fig. [3j because the coloured 
domains are clearly much larger than the spacing between 
(black) pinned sites. In other words, each correlated re- 
gion in these images contain very many pinned spins. 



These observations will be quantified below in Sec. IV C 



The qualitative description of these images is strongly 
reminiscent of observations made in dynamic heterogene- 
ity studies |10j . except that time has now been replaced 
by the fraction of pinned spins. The images in Fig. [^sug- 
gest that a similar behaviour is found in both square and 
triangular models but that the specific features of the 
models will be reflected in the form of the correlation 
functions. For example, the SPM is characterized by 
strongly anisotropic correlations, while correlations ap- 
pear more isotropic for the TPM, although they do have 
an underlying fractal structure. These observations once 
again echo previous studies of the dynamic heterogeneity 
in these models [16] 



B. SPM: 'Bulk' scaling at / = 

To analyse length scales and their scaling in the SPM, 
it is useful to start by considering static correlations for 
the 'bulk' at / = 0. As discussed in Sec. [II] the first 
non-trivial correlations involve four spins arranged at the 
edges of a rectangle of size x x y. It is clear from Eq. ^ 
that lines xy = const, are contours of the static four- 
point function g s {x, y). That is, the four-spin correlations 
are strongly anisotropic, which leads to unusual scaling 
behaviour at low temperature (c — > 0) . For example, one 
may measure correlations at a fixed finite angle to the 
lattice axes (with ^ 0, 7r/2, etc), in which case 



g s (x,y) ~ G so (ry/c,0), 



(11) 



with x — r cos and y = r sin 0, as usual. (Explic- 
itly G so (u,0) — e~ u sm2e .) Since the scaling variable 
is u = r-\fc, the correlation length away from the lat- 
tice axes scales as l/y/c. (Throughout this section we 
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use the symbol G for scaling functions, with the approxi- 
mate equalities valid on taking c — > with the arguments 
of G held constant.) 

However, a larger static correlation length in this sys- 
tem is revealed by measuring g a (x,y) along the axes of 
the square lattice. For fixed y (of order unity) and vary- 
ing the temperature, one gets 



s{x,y) - G sa {xc, y) 



indicating a correlation length £ ~ c 
the lattice axes. (Explicitly G sa (u,y) 



(12) 

measured along 



_ P -2y« 



At low temperatures, one may also show that the cir- 
cular average of g s (x,y) is dominated by contributions 
from near the axis and so it also decays on a lengthscale 
c" 1 , as g s (r) ~ r~ 1 G sr (rc), which may be rewritten 



as 



g s (r) ~ cG sc (rc). 



(13) 



That is, the effect of the circular average is to pick up 
the longest of the two lengthscales that appear in g s (r). 



C. SPM: Real-space scaling at / = /* 

In Fig. |4j we show the behaviour of g^j and Xij for 
/ = /* and a representative temperature /3 = 3. We 
compare these correlations with the behaviour of g^ij (t) 
measured at / = and t = t, for the same tempera- 
ture. In all cases, the correlations seem to operate over 
a similar lengthscale (the same linear scale is used for 
all panels). The correlation functions are all strongly 
anisotropic, although we observe slightly different angu- 
lar dependences in each case. 

The dominance of a single lengthscale in this problem 
may be seen from Fig. [5] where we show circular averages 
of these correlation functions, plotted as a function of 
the scaling variable rc. It was shown in Ref. |16j that 
dynamical four-point correlations at / — and t = r 
collapse as a function of this scaling variable. Here we 
show that the same behaviour holds for and g^ij 
evaluated at /*. 

To describe the scaling behaviour of correlation func- 
tions near /* in the SPM, we make an ansatz for the 
circular averaged three-point function X{r): 



X(r)\ r ^cf*G x (rc). 



(14) 



The choice of scaling variable rc indicates that the dom- 
inant correlation length in the system scales as £ ~ cr 1 . 
To understand the prefactor c/*, note that the ratio 
Xij/ f quantifies the effect of freezing spin j on the au- 
tocorrelation function at site i. Thus, if each frozen spin 
has an O(l) effect then one would expect X(r) fG(r), 
which explains the presence of /* as a prefactor in 
Eq. (fl4| . The extra prefactor of c in Eq. ( 14 ) has two 



possible interpretations, which are hard to discriminate 
on the basis of our numerical results. First, if the correla- 
tion function X{x,y) is largest near the lattice axis, and 
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FIG. 4: Real-space correlations in the SPM. (a) The correla- 
tion function Xij evaluated at / = /*, in the long-time limit. 
The co-ordinates x and y give the position of spin j relative 
to spin i. (b) The four-point susceptibility g^^j at / = /* 
and long times, (c) The four-point susceptibility at / = 0, 
evaluated at the bulk relaxation time r. The scales are the 
same in all plots, indicating that all correlations operate over 
similar length scales. 



if these on-axis correlations are 0(1) and dominate the 
circular average then one arrives at X(r)/f rs cGx{r), 
as in the case of the static function g s {r). However, a 
second explanation could be the presence of off-axis cor- 
relations of strength 0(c), which would lead to the same 
prefactor c in Eq. ( 14 ) . 



Our ansatz for the low-temperature scaling of of g^ is 

g<(r)\ f . «/*(l-/*)G 4 (rc). (15) 

Again, the scaling variable rc indicates that the corre- 
lation length scales as £ ~ c _1 . In Fig. [5j we plot 
<74/[/*(l — /*)] as a function of the scaled variable rc, 
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FIG. 5: Circular averaged correlation functions in the SPM 
evaluated at / = /* and plotted to illustrate their scaling 
with temperature, (a) Circular average of Xij. (b) Circular 
average of g^jj . In both cases the relevant length scales scale 
as £ ~ c _1 . At the lower temperatures we show data points 
are only for those r and ft where numerical uncertainties are 
smaller than or comparable to symbol sizes. 



a procedure which nicely collapses our data. The pref- 
actor (1 — /*) is irrelevant for the purpose of scaling in 
the low temperature limit. However, it is natural from 
a physical point of view because g^n oc (1 — /), and we 
do find that it improves the data collapse in the stud- 
ied range of temperatures. The physical interpretation 
of the prefactor /* in Eq. ( |15|) i s not immediately clear. 
We note that Eqs. (14 1 and (1151) are together consistent 



with Eq. ( 10 1, which holds if the correlations of the are 
directly attributable to individual frozen spins fj. More 
support for this can be obtained through a direct numer- 
ical evaluation of the right hand side of Eq. ( 10 ) which 
has the same dependence on ij as g^ij, but is smaller by 
factor close to 5, independently of the temperature. As a 
result, xa & n d its bound A% 4 scale in the same way, but 
differ by a prefactor. 

To conclude, we have shown robust evidence that g& %j 
and Xij in the SPM are both controlled by the same 
lengthscale £ ~ (V c )- The scaling of the prefactors in 
these correlations is less clear, but Eqs. ( 14 1 and ( 15 ) 
are consistent with our numerical data. Assuming that 
these results do hold, we arrive at the following scaling 
behaviours for the susceptibilities: 



df 



Ax- 



41/* 



(16) 



We recall that dynamical correlations (74 (r, t)|/=o a t 
t = t are controlled by the same lengthscale £ ~ (1/c), 
but that Xi(t — T )l/=o ~ c due to the strong 
anisotropy of the correlation function. Combined with 
the super- Arrhenius growth of the relaxation time shown 
in Fig. [if, this difference in the scaling of \4 emphasizes 
that the relaxation near /* is qualitatively different from 
bulk relaxation at / = 0, even if the same lengthscale ap- 
pears in both cases. In particular, the susceptibility X4 at 
/* grows more quickly on cooling than the bulk X4j con- 
sistent with the observation that relaxation is slower and 
more co-operative. In this respect increasing /* is sim- 
ilar to reducing the temperature. However, in contrast 
to decreasing T, there is no evidence for an increasing 
length scale as / is increased. 



D. TPM: Real-space scaling at / = /* 

For the triangular plaquette model, the isotropic im- 
ages shown in Fig. [3] lead us to compute circularly aver- 
aged correlation functions directly. The results are pre- 
sented in Fig. [6j showing scaling with temperature. We 
find that our data are most consistent with 



and 



X(r)\ r « crG XT (rc), 



94(r)\ r «/*(!- f*)G 4T (rc). 



(17) 



(18) 



The most striking feature of Eqs. ( 17 1 and ( 18 ) is that 



the scaling variable rc is the same as that found in the 
SPM. This suggests that correlations at /* extend over 
a length which scales as £ ~ 1/c. This is surprising be- 
cause the point-to-set length for the TPM at / = does 
not scale as cr 1 , nor does the dynamic correlation length. 
For dynamical correlations at / = and t = t then it 
is known [T7] that ^(r, t)|/=o ~ G , 4T(rc 1 / df ), such that 
both static and dynamic lengths scale as £(/ = 0) ~ 
c~ Y l di ~ c~ - 631 . In Fig. [g]:, we show that these two 
scaling forms can be clearly differentiated over the tem- 
perature range shown, since g^{r)\(—Q does not collapse 
as a function of rc, as expected. It is therefore clear 
that a new lengthscale £ ~ c _1 appears in the TPM near 
/ = /*, which is longer than any static or dynamic cor- 
relation length found previously for the bulk at / = 0. 

As in the SPM, the physical interpretation of the scal- 
ing prefactors in Eqs. (17) and (18) is not clear. The 



scaling laws we have proposed indicate that the bound 
A%4 scales in the same way as x 4 , although A%4 is sig- 
nificantly smaller than \4 in the TPM. We also note that 



while the scaling forms in Eqs. ( 17 1 and ( 18 ) in the TPM 



/* 



are the same as in Eqs. ( 14) and ( 15 ) in the SPM, we do 
not see any a priori reason for this result. In particular, 
the spatial structure of the correlations are quite different 
in both cases. 

While the appearance of a new lengthscale near /* 
makes the TPM different from the SPM, we emphasize 
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FIG. 6: Circular-averaged correlations in the TPM. (a) Cir- 
cular average of Xij at /*, showing a length scale £ ~ c" 1 . 
(b) Circular average of at /*, showing a similar length 



scale £ 



(c) Circular average of gi,ij(t) evaluated at 



f — and t — t. This correlation function clearly has a dif- 
ferent scaling to those evaluated at /*: it was shown in [IT] 
that £ ~ c _1,/df in this case. 



that the relaxation mechanism changes qualitatively near 
/* in both models. In the SPM, this appears as a larger 
relaxation time and a larger susceptibility without any 
increase in the length; in the TPM the lengthscale, time 
scale and susceptibilities are all different from their values 
at / = 0. 



V. SQUARE PLAQUETTE MODEL: 
AND STRONG PINNING 



WEAK 



We have discussed the scaling of length and time scales 
at /*, as temperature is reduced. In this section, we 



(a) 





4J(= k) 



FIG. 7: Two diagrams representing contributions to Xij in 
the SPM, at order / 3 . The white circle represents spin i, the 
grey circle spin j, and the positions of spins klm are to be 
summed, subject to the constraint that j is equal to one of 
klm. (a) Leading contribution to X(x,y), assuming x,y > 
0. (b) Leading contribution to X(x,0): summing over the 
positions of spins / and m corresponds to a sum over the 
dimension y' . 



consider how /* depends on temperature. We focus on 
the SPM for which analytic calculations provide useful 
insight. 



A. Small-/ limit 



We concentrate on the behaviour of the correlation 
function X^. We define a parameter fj, by / = (1 + e M ) _1 
so that the limit of small / is equivalent to a limit of small 
e~ M . The correlation function has a series expansion 
in powers of e _/i given by 



Xi 



1 



k<l 



?/(!-/), (19) 



where Zi = (l+c~'') y and (• • • is an average in a sys- 
tem where spins (k, I, . . . ) are pinned and all other spins 
are mobile (unpinned). For details, see Appendix |A| par- 
ticularly Eq. (lA4|. The factor of q in Eq. (fl9|) must be 



obtained by a separate series expansion over /. 

Assuming i ^ j, symmetries of the SPM imply that 
the first non-zero term in the expansion is at third order, 



-3/i 



k<l<m 



j I klm ■ 



(20) 



The factor fj means that the average is zero unless j is 
equal to one of k, I, or m. As shown in Fig. [7[ the corre- 
lations may be calculated in a diagrammatic expansion. 
Spin i is shown as a white circle and has a fixed position. 
Spin j is shown as a grey circle: its position is fixed, and 
we also have the constraint that one of the frozen spins 
k, I, or m coincides with j. Equation (191 shows that 



we must sum over the positions of the remaining frozen 
spins: these spins are shown as black circles. To eval uate 
the contribution of each diagram to Xu , we use Eq. ( A5 ) 
which shows how to evaluate expectation values in the 
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presence of a fixed set of frozen spins. For the SPM, 
(a-ifj)kim is non-zero only if spins iklm lie on the four 
vertices o f a rectangle. If the rectangle is of size x x y 
then Eq. yield s (aifj)klm = tanh(^/2) 2 ^. (The de- 
nominator in (A5| has a trivial value 2 -3 in this case, 



(a) 



because all configurations of the frozen spins are equally 
likely and they all have equivalent effects on spin i.) 

The leading order behaviour of Ay stems from two 
distinct cases, as shown in Fig. [7j If spins i and j are in 
the same row of the square lattice, with spacing x, then 
we fix k — j and we sum over all sites I and m such that 
the sites iklm fo rm a , rectangle. For a rectangle of size 
x x y' then Eq. (Xi) yields {cnfj) Hn = tanh(/3/2) 2 ^'. 
Summing over the positions of spins / and m, one obtains 
a geometric series and the result is, for x ^ 0: 



A(x,0) = / 3 



tanh(/3/2) 2 N 
1 - tanh(/3/2) 2 N 



0(f 



(21) 



where we have defined X(x, y) — Xy, evaluated for spins 
i and j that are separated by a vector (x,y). A similar 
analysis applies if spins i and j are in the same column 
of the lattice, yielding A(0, y) — X(y, 0). 

However, if spins i and j are in different rows and 
columns then we fix I — j and we sum over sites k and 
m such that sites iklm still form a rectangle. There 
is only one choice for k and m in this case, as shown 
in Fig. [7)3. If the vector from site i to site j is (x, y) 
then the resulting rectangle is x x y in size and, again, 
(aifj)kim = tanh(/3/2) 2xy so that, for x,y ^0, 

X(x, y) = f tanh03/2) 2 ^l + 0(/ 4 ). (22) 

Collecting all these results and summing over the vol- 
ume, Eq. ([I]), one finally obtains the leading order be- 
haviour of the three-point susceptibility, namely 



Of 



12 A c (T)f + 0(f), 



(23) 



where 



tanh(/V2) 2 - 
c[ ' 4-f 1 - tanh(/3/2) 2 * 



-[/3 + 0(l)]. (24) 



The final approximate equality holds for small c (i.e., at 
low temperature) and follows because YlkQ — /(l — 
{1-S) k ) « (l/<5)[log(l/<5) + 0(l)] as 5^ 0. InFig.[§i we 
show that the result in Eq. ( 23 ) holds very well for small 
values of /a 



!A C , but breaks down for larger /, where 
higher order terms in the expansion also contribute, as 
discussed below. 

The physical interpretation of this small-/ result is 
that adding one extra frozen spin j affects the autocor- 
relation a,i over an area of linear size £ ~ c , but the 
anisotropy of correlations imply that the correlation vol- 
ume is v ~ A c ~ (— In c) / c, much smaller than the naive 
assumption v ~ £ 2 ~ 1/c 2 . The strength of the response 
on adding the frozen spin is small (proportional to / 2 ) 




(b) 
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FIG. 8: Plots of q and 9g/9/ for the SPM, as a function 
of the scaling variable jyfA~ c . (a) Plot of the overlap q in- 
dicating that /* ~ Ac 1 ^ 2 . (b) Slight change in data. The 
susceptibility dq/df i s m aximal near /*. The straight line is 
the theoretical result (231 for the small-/ limit. This approx- 

lA~c <g 1 and does not capture 



imation applies only when f\ 
the peak in dq/df near /*. 



but the lengthscale controlling Xy is large and indepen- 
dent of / as / — > 0, and the sum in Eq. (|8| is dominated 
by correlations close to the lattice axes. 



Using Eq. (23), it is then easy to integrate dq/df to 



obtain the low-/ behaviour of the overlap, 
q = 4A c f + 0(f i ). 



(25) 



It is interesting to remark that A c = Yl X y(>o) 9s( x i v) 2 > 
which highlights the fact that in the limit where the ran- 
domly frozen spins are dilute, the static overlap and sus- 
ceptibilities simply capture the most trivial behaviour of 
the bulk system. In a supercooled liquid where two-body 
correlations do not vanish, one would expect the overlap 
to be proportional to / for small-/, with a prefactor di- 



rectly given by the pair correlation function, see Eqs. ( 19 ) 
and (A8). It is only by going beyond the leading order 



in / one reveals the relevant higher-order correlations re- 
sponsible for amorphous order [ITJ [23] . 



B. SPM: Behaviour near / = /* 

To analyze the behaviour near / = /*, we now con- 
sider higher-order terms in the expansion over /. In 
Fig. [9j we show two further contributions to Ay , in the 
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FIG. 9: Two representative diagrams that contribute to Xij 
in the SPM, at order f 5 (left) and f (right). The dimensions 
x and y are fixed by the positions of spins i and j while the 
primed dimensions are to be summed over. These sums yield 
factors of A c (left) and A\ (right) in the contributions of these 
diagrams to Xij. 



case where i and j are not in the same row or column. 
These contributions appear at order e _5M and e~ 7AI in 



Eq. 



(19 1, respectively. For these diagrams, Eq. 

= tanh(/3/2) 2 - 4 where A is the 



(A5) 



yields (ai/j)fc;... = tanh(/3/2) 2 ^ where A is the total 
area of the rectangular regions enclosed by the solid 
lines in the diagrams. (The denominator in Eq. (A5) 
is again trivial for these diagrams, due to symmetries 
of the SPM.) Summing over the positions of the frozen 
spins yields a factor of A c for each rectangle, except for 
the rectangle whose location is fixed by the positions of I 
and j. These two diagrams therefore contribute to as 
f 3 (f 2 A c )t&nh((3/2) 2x y and f(f 2 A c ) 2 tanh(/3/2) 2x », re- 
spectively. Of course, there are various other contribu- 
tions at these orders that scale in the same way. Con- 
structing higher order diagrams similar to those in Fig. |9j 
one may identify a series of positive terms which are all 
proportional to powers of f 2 A c . For small c, we expect 
these to be the largest terms at each order. These con- 
siderations clearly motivate the use of f\fA~ c as a scaling 
variable for this expansion. Of course, if f 2 A c is not 
a small number, then the leading terms in the small-/ 
expansion do not give a good approximation to the cor- 
relation function. 

In Fig. [8] we plot q and dq/df as a function of the 
scaling variable f-J A C (T). We find that the crossover /*, 
which corresponds to the maximum observed for dq/df, 

— 1/2 

does indeed scale as A c , so that finally 



f*(T) 



In c 



Texp 



1 

2T 



(26) 



This result shows that the 'localization' crossover at /* 
occurs along a line /* (T) in the (/, T) phase diagram, 
with /* -> as T -> 0. 

Also, the numerical data indicate that results from the 



small-/ limit such as Eq. (23) are applicable only when 



'Ac <C 1 and break down for / ~ /*. As suggested 
above, this strongly suggests that the maxima in dq/df 
and %4 have their origin in non-trivial "many-body" ef- 



fects that are not captured by the low-order expansion 
about the dilute limit. 

Physically, the interpretation is that the relaxation 
mechanism near to /* is qualitatively different from the 
bulk relaxation at / = 0. Perturbation theory in / is 
not sufficient to capture this new mechanism: a non- 
perturbative approach would be necessary to make fur- 
ther analytic progress, presumably by summing infinite 
subsets of diagrams in the expansion of correlation func- 
tions. In this model, /* — > when T — »• 0, indicating 
that the regime where perturbation theory is valid be- 
comes vanishingly small near the glass transition of the 
model (which takes place at T = 0). 



VI. DISCUSSION 

We have studied the effect of random pinning in the 
context of two finite dimensional spin models with pla- 
quette interactions. By increasing the fraction / of 
pinned spins at a fixed low temperature, we have discov- 
ered the existence of a temperature-dependent crossover 
between bulk- like relaxation at small / < f*(T), and 
a nearly localized amorphous state at large / > f*(T). 
The study of static correlation functions and susceptibil- 
ities in the presence of random pinning directly reveals 
the existence of growing amorphous order on cooling. 
This growth appears through the large lengthscales that 
can be measured by and g& t ij, and by the decreas- 
ing values f*{T) required to keep the system localized 
in a single state at lower temperature. Moreover, these 
measurements do not require a priori knowledge of the 
specific type of order that sets in at low temperature. 
These results therefore demonstrate that the main ob- 
jective underlying the measurement of point-to-set cor- 
relation functions is fulfilled in plaquette models. 

The lengthscales that we measure show scaling be- 
haviour at low temperatures. We have emphasized that 
while low-order terms in the expansion over / are re- 
lated to static correlation functions of the bulk system at 
/ = 0, the length and time scales that we observe near 
/ = /* are related to nonperturbative effects, and an- 
alytic calculations of lengthscales in that regime would 
seem to require a new approach beyond those given here. 

Similarly, we have shown that length and time scales 
near /* in the plaquette models are not related in the 
same way as they are at / = 0. Increasing / in the 
SPM, we found an increasing time scale and a growing 
susceptibility, but without any increasing length scale. In 
the TPM, the length scale increases as / is increased but 
the change in time scale is very mild, in contrast to the 
strong dependence observed at / = for this model. 

We have also emphasized that although static and dy- 
namic correlations are strongly enhanced near /* as com- 
pared to the bulk, length and time scales are finite at /* , 
after taking the thermodynamic limit at any non-zero 
temperature. This implies that random pinning does not 
induce any kind of phase transition in plaquette models. 
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It is perhaps unsurprising that these models do not ex- 
hibit an ideal glass transition in the (/, T) phase diagram, 
since no transition occurs at finite temperature in the 
bulk at / = either. However, the sharp crossovers we 
have revealed in plaquette models represent non-trivial 
new results, because they have no counterparts in the 
bulk systems at / = 0. 

Given that neither mode-coupling theory nor RFOT 
theory represent accurate descriptions of plaquette mod- 
els in the bulk, we do not expect these approaches to ac- 
count for the effect of random pinning either. Thus, we 
argue that the results obtained within plaquette models 
provide a useful alternative reference point for interpret- 
ing simulation data for more realistic models of super- 
cooled liquids. In particular, the absence of any phase 
transition at /* indicates that such transitions may not 
be generic in glassy systems with pinning. 

In this respect, it is instructive to compare the re- 
sults we find here with the predictions of RFOT the- 
ory [T2]. A central quantity in this theory is the config- 
urational entropy, which measures the diversity of long- 
lived metastable states. Assuming that RFOT applies in 
supercooled liquids, results for model systems and renor- 
malisation group calculations |12j indicate that the con- 
figurational entropy density s c is well represented by 



s c (T,/)«s c (T,0)-/Y(T). 



(27) 



In 3 dimensions and above, this leads to a phase transi- 
tion at f*(T) w Y(T)/s c (T,0). In addition, RFOT pre- 
dicts that a lengthscale £pts grows as s c (f, T) approaches 
zero, and that the relaxation mechanism at finite / in- 
volves cooperative rearrangements over the lengthscale 
£pts, much as in the bulk. We emphasise that s c is re- 
lated to metastable states and may not be obtained from 
the statistics of minima on the system's energy landscape. 
A precise definition of s c in finite-dimensional systems is 
slightly problematic since all metastable states have fi- 
nite lifetime in that case. However, one may follow the 
procedure of Ref. [5] , as long as the timescale associated 
with structural relaxation is well-separated from all mi- 
croscopic timescales. 

Turning to the plaquette models, the geometric con- 
struction of Ref. [H) indicates that s c decreases as / in- 
creases, just as in RFOT theory. However, it is clear from 
Fig. [3] that even for / > /* there are sets of spins that 
may rearrange cooperatively, which ensures that s c does 
not vanish at /*. Thus, while s c presumably decreases 
sharply near /*, it does not drop to zero as predicted 
by Eq. (27). For supercooled liquids, RFOT predicts in- 
stead localization in a single state for / > /*, so that if 



the analysis of Sec. IV were repeated for those systems 
then regions where cooperative motion is possible should 
be forbidden for / > /*, and the light-coloured regions 
shown in Figs. [3] would be completely absent. In plaque- 
tte models, the existence of such regions restores a finite 
configurational entropy density above /* and the pro- 
posed phase transition is avoided. It is unclear whether 
such strong spatial fluctuations can be present in super- 



cooled liquids, and whether they are properly captured 
by renormalization group treatments [12 . This remains 
an area of ongoing research 

Thus, the plaquette models illustrate that even if the 
detailed predictions of RFOT do not apply, systems 
where configurational entropy decreases on pinning can 
be generically expected to exhibit increased cooperativ- 
ity on increasing /, accompanied by growing timescales 
as well as growing lengths and/or growing static suscep- 
tibilities. One may also expect crossovers that sharpen 
and move to small / on cooling. To this extent, the 
plaquette models are broadly consistent with published 
numerical results for particle models [TTJ l2"D - f2"3l [33] , al- 
though the /-dependence of the relaxation time appears 
much weaker in the plaquette models in comparison with 
model liquids. 

On the other hand, the most striking prediction of 
RFOT theory is the presence of a phase transition at fi- 
nite /*. This transition remains to be found numerically, 
as the behaviour for large values of / and low temper- 
atures has not been investigated in much detail so far. 
In the absence of such a transition, we argue that the 
detailed RFOT scaling predictions for length and time 
scales must be tested directly in order to substantiate the 
theory and distinguish it from a more general picture of 
increasing cooperativity in the presence of pinning. Test- 
ing these predictions remains however a very challenging 
task, especially in the absence of direct measurements of 
the configurational entropy. We therefore conclude that 
while pinning particles is an interesting new method of 
measuring amorphous order and its growth upon cooling, 
it does not necessarily resolve the central problem of how 
to test the fundamental assumptions of RFOT theory by 
practical measurements. Nevertheless, we hope future 
studies will investigate further the effect of random pin- 
ning in supercooled liquids, especially in the relatively 
unexplored regime of strong pinning. 
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Appendix A: Correlations in systems with random 
pinning 

In this appendix, we discuss some general results for 
spin systems in the presence of pinning. 



1. Ensemble dependence of x*(f) 

In the systems considered here, the set of frozen spins / 
remains constant as the dynamics proceeds. This results 



13 



in ensemble-dependent susceptibilities 34-36], such as 

X4(t)=X4F(t) + A X 4{t), (Al) 

where X4f(*) = [(1 — /) V] — 1 (i) 2 ) / is evaluated in a 
'restricted' ensemble with a fixed number of frozen spins. 
Analysis of such ensemble-dependence can be useful for 
understanding how the time-independent variables (the 
frozen spins in this case) influence the time-dependent 
ones. 

The difference term Ax4,(t) may be derived as in 
Ref. [3BJ or equivalently following Ref. [33]. We write Q — 
(Q) = (Q - (Q)f) + ((Q) f-(Q)) where (•••)/ is an aver- 
age with fixed Nf, as above. Substituing into Eq. ([H]), we 
note that if the restricted ensemble has Nf ~ fV frozen 
spins then (Q) f - {Q) « V~ 1 (d{Q) /df)(Nf - fV) which 
gives (5Q{t) 2 )^ = {5Q{tf) f + V- 2 (0(Q) /df) 2 ((6N f ) 2 ) 
(the equality is exact in the limit of large system size 
V, since the fluctuations of Nf are small in that case). 
Noting that ((5Nf) 2 ) = Vf(l - f) and d(Q)/df = 
V[(l - f)dC(t)/df - C(t)} then the result ^ follows. 

Since XiF(t) and Ax4(t) are both non-negative then 
Xi{t) > Axi{t). If this bound is saturated, this means 
that the correlations between the a, are directly at- 
tributable to the influence of individual fj. In partic- 
ular, a sufficient condition for Xiify — A X 4,(t) is that 
for a fixed choice / of frozen spins, the autocorrela- 
tions di are all independent and respond linearly to 
the fj, so that (ai(i))/ = {ai(t)) + J2 k fkU lk (t) where 
Uik(t) = Xik(t)/[f(l — /)] is assumed independent of /. 
Physically, this condition means that the site-to-site fluc- 
tuations of the a, depend only on the frozen spins, and 
the effect of each frozen spin is independent. 

In addition, independence of the (a«) / means that 
(ai(ij)f = (cii) f(cij} f. Following Ref. [37] we use an over- 
bar to indicate the average over /, so the definition of the 

-, ; — 2 



(did 



ilf 



Hence. 



four-point function is g^i 
94,ij = I2kk>[fkfk> - f 2 ]U lk U jk ,, and si nce th e frozen 
spins are chosen independently one has fkfk' ~ f 2 — 
f(l — f)5kk'- The resulting expression for g^j is given 
in Eq. (10) of the main text, where g^^j appears as a 



convolution of with itself. 

We note that for Ising spin variables Sj, it is not pos- 
sible to satisfy (didj)f — (di}(dj)f in the case i = j. 
However, assuming that Xi i s dominated by collective 
behaviour and not single-site fluctuations, one still ex- 
pects Eq. ( 10 ) to hold as an approximate equality if the 
di are primarily determined by a linear response to the 



2. Long time limit 



As discussed in Sec. |III[ long-time limits of correla- 
tion functions in systems with pinned spins are static 
(thermodynamic) quantities and can be calculated within 



equilibrium statistical mechanics. The analysis is similar 
to that of Krakoviack [37] for particle systems. 

We write spin configurations as s = (s±, . . . , sy) and 
choices of frozen spins as / = {fi,---,fv)- The dis- 
tribution of the initial (reference) configuration s A is 



Pi(s A ) 



1 jZ\ where E(s ) is the energy of con- 



figuration s A and the partition function Z\ enforces nor- 
malisation. The fi are all independent with (/,•) = / and 
so their distribution is Pf(f) = e^^^i** /Zf where [i is 
defined through f — (1 + c' 1 ) -1 and Zf is a normalisa- 
tion constant. We denote the number of frozen spins by 
Af = Ei fu noting that Nf + N m = V. 

In the long-time limit, and for fixed s A and /, the final 
configuration s B has probability distribution 



P 2 (s B |/,s A ) 



^b(/,* a ) 



-^HKi-fM^}, 

(A2) 

where Zs{f, s A ) is a normalisation factor, defined so that 
Es B -f > 2('S B |/, s A ) = 1. Thus, in the long-time limit, av- 
erages (• • • ) are taken with respect to the distribution 



P(f,s A ,s") = F f (/)A(s A )P 2 (^|/, S A ). (A3) 

It may be shown that this distribution is invariant under 
s A <-> s B . In particular, this means that the marginal 
distribution of s B is equal to the Boltzmann distribution 
Pi(s B ). Physically, this means that structural averages 
and correlation functions arc unaffected by the pinning. 
Long-time correlations between the fi and the di may be 
calcu lated as averages with respect to the distribution 
(A3), identifying di = (1 - fi)s A sf. 



3. Small-/ limit 

Correlation functions between the di and the fi in the 
long-time limit may be calculated in an expansion about 
/ = 0. The idea is simply to collect together configu- 
rations / with exactly 0, 1 ,2,. . . frozen spins. Formally, 
one expands Pf(f) in Eq. (A3) over e~ M , so that for any 
observable Y, 



(Y) 



1 

Zf 



(A4) 

where (Y)jki... is an average over configurations A and B, 
given that spins jkl . . . are frozen and all other spins are 
unfrozen. (The average (Y)q is taken without any frozen 
spins, so that A and B are independent configurations 
from the Boltzmann distribution and averages factorise 
as (h(s A )f 2 ( S B )) = (/ 1 (^))o(/ 2 ( S s ))o.) 

To make progress, the key st ep is to write the individ- 
ual expectation values in Eq. (A4) as correlations with 
respect to the distribution Pi(s A )Pi(s B ), in which case 
configurations A and B are independently chosen config- 
urations from thermal equilibrium (at / = 0). We write 
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(• ■ • )o for averages with respect to this distribution. Us- 
ing Eq. (A3), the general result is 



(Y) k 1 k 2 ...k p 



Y\ 



2-"IK = i(i + «) ' 

Z B (f,s A ) 



(A5) 



where the notation Y\f indicates that any /-dependence 
of Y\f has been accounted for by substituting the specific 
set of frozen spins k\ . . . k p . For example if Y = ciifj as 

in the main text then Y\ f a i^2r=i k r since Y = 

unless fj = 1. To obtain (A5) we used 5 SS > = |(1 + ss') 
for Ising spins s, s' , and we also have that 



M/,o=2-*(n(i+<*o 



where the average is taken over the spins s, i.e., 
Es( ... )e ^)/T /Zi . 



(A6) 



For example, Eq. (A4| can be used to calculate q = 
(ttj), leading to 



Q = 




-V( 



-2fj, 



), (A7) 



where we used Z{ = (1 + e"^) y and we assumed that 
averages are translationally invariant, (s,) = (sj)o etc. 

If the system has inversion symmetry so that (sj)o = 
then the denominator in Eq. ( |A7[ ) is trivial and the av- 
erage may be evaluated directly. Physically, the denomi- 
nator accounts for the fact that the autocorrelation (dj) f 



depends on the state of s, in the reference configuration 
s A , and the different values of sf may not be equally 
likely. In the symmetric case, both values of sf are 
equally likely so the denominator has a trivial value. 

To make progress with the general case, note that Si = 
±1 which means that that for any function F(s) then 
F(s) = |(1 + s)/(l) + |(1 - s)f(-l) and hence (for any 
x), 1/(1 + xsi) = (1 — xsi)/(l — x 2 ). The result is that 



MS 



'El 



(SsiSs 



3/0 



s i)o 



0(f 



(A8) 



One may similarly show that 



, (SsiSs 



3/0 



1 



s i)o 



o(/ 2 ), 



(A9) 



and summing over j and integrating with respect to / 
yields Eq. ([A8), via Eq. 



Physically, Xij/ f is the change in the autocorrelation 
(ttj) if one restricts to an ensemble where fj — 1. If the 
spin syst em has two-point correlations of range £ then 
Eq. ( A9 1 shows that pinned spins influence mobile spins 



over a range of at least £, resulting in an O(f) contri- 
bution to q. Of course, if two-spin correlation functions 
dominate the physics then the pinning procedure is re- 
dundant since the correlations may already be measured 
through the spin-spin correlation function. For the spin 
models we consider in this paper, multi-spin correlations 
dominate the physics, and we find that X$j is a useful 
way to reveal the relevant correlation lengths without re- 
quiring explicit measurement of multi-spin correlations. 
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